Airway disease decreases the therapeutic potential of epithelial stem cells

Backgorund Tissue-engineered tracheal grafts (TETG) can be recellularized by the host or pre-seeded with host-derived cells. However, the impact of airway disease on the recellularization process is unknown. Methods In this study, we determined if airway disease alters the regenerative potential of the human tracheobronchial epithelium (hTBE) obtained by brushing the tracheal mucosa during clinically-indicated bronchoscopy from 48 pediatric and six adult patients. Results Our findings revealed that basal cell recovery and frequency did not vary by age or region. At passage 1, all samples produced enough cells to cellularize a 3.5 by 0.5 cm2 graft scaffold at low cell density (~ 7000 cells/cm2), and 43.75% could cellularize a scaffold at high cell density (~ 100,000 cells/cm2). At passage 2, all samples produced the number of cells required for both recellularization models. Further evaluation revealed that six pediatric samples (11%) and three (50%) adult samples contained basal cells with a squamous basal phenotype. These cells did not form a polarized epithelium or produce differentiated secretory or ciliated cells. In the pediatric population, the squamous basal cell phenotype was associated with degree of prematurity (< 28 weeks, 64% vs. 13%, p = 0.02), significant pulmonary history (83% vs. 34%, p = 0.02), specifically with bronchopulmonary dysplasia (67% vs. 19%, p = 0.01), and patients who underwent previous tracheostomy (67% vs. 23%, p = 0.03). Conclusions In summary, screening high-risk pediatric or adult population based on clinical risk factors and laboratory findings could define appropriate candidates for airway reconstruction with tracheal scaffolds. Level of evidence. Level III Cohort study. Graphical Abstract


Background
Laryngotracheal stenosis (LTS), defined as soft tissue narrowing of the proximal airway, currently has an estimated incidence of 0.6-2% [1][2][3].The causes of LTS can be multifactorial in nature.In the pediatric and adult populations, acquired stenosis is most commonly secondary to prolonged intubation, tracheostomy, systemic disease, trauma, infection, or idiopathic causes [4,5].Surgical management of LTS can range from tracheostomy to endoscopic and open approaches [1].Despite advancements in endoscopic and open techniques, there are significant limits to current clinical approaches.
Long-segment tracheal defects are defined as > 50% in adults and > 30% in children [6].These defects cannot be repaired primarily because they require replacement tissue that is not currently available.A myriad of

Graphical Abstract
bioengineering approaches have created scaffolds that recapitulate the biological and mechanical properties of the native trachea [7,8].The success of these scaffolds is dependent in part on rapid epithelialization, which restores mucociliary clearance and prevents infection and inflammation [9].However, the optimal method for reepithelialization is unknown.
We demonstrated that host-derived basal cells regenerate a surface airway epithelium (SAE) on a partially decellularized tracheal scaffold [10].This process requires 7-14 days and could be accelerated by preseeding the scaffold with patient-derived cells.Ex vivo methods of cell seeding include the use of induced pluripotent stem cell-derived basal cells and primary airway stem cells (basal cells) [11].Although both cell types can be expanded using current cell culture technology, the quality of the cell therapy product is likely to be dependent on the initial cell inoculum.
Patients with long-segment tracheal defects can suffer from chronic airway disease, a situation that may compromise the regenerative potential of the host's basal cells.In particular, stenotic regions of the airway exhibit abnormal epithelial morphology, including basal cell dysplasia [12,13].We reported that the lifespan of the tracheobronchial basal cell is ~ 40 population doublings [14], which is consistent with the Hayflick limit [15].Stem cell frequency decreases during wound healing and in vitro expansion, resulting in reduced regenerative potential [16,17].These decrements are associated with a decrease in the number of clone-forming stem cells, accumulation of cells with shortened telomeres, and suboptimal production of differentiated secretory and ciliated cells [18].Thus, a critical question is whether tracheobronchial stem cells from patients with airway disease have sufficient regenerative potential to support the cellularization of a graft and the production of a functional epithelium.We addressed this issue by assessing the phenotype and function of basal cells that were recovered from the airways of normal donors and those with chronic lung disease.

Approvals
This study was approved by the institutional review boards of two tertiary-level pediatric and adult hospitals (Nationwide Children's Hospital IRB STUDY00000847, Ohio State University 2021N0027).Pediatric and adult patients undergoing scheduled direct laryngoscopy and bronchoscopy (DLB) were voluntarily recruited.Demographics of the patients recruited were recorded.

Tracheal airway epithelial cell isolation
Our methods are similar to those previously described [16,18].Briefly, endotracheal brush biopsies were collected using a cytology brush (Medical Packing Corporation, Camarillo, CA).Under direct visualization, the cytology brush was advanced and spun in a focal region of the trachea five times for approximately 10 s.

Tracheal airway epithelial culture timeline
Cells were cultured using the modified conditional reprogramming culture (mCRC) method [19].Accordingly, the cells were co-cultured on a feeder layer of irradiated NIH-3T3 fibroblasts (ATCC #CRL-1658) with F-medium (F med ) supplemented with Rho-kinase inhibitor, Y-27632 (10 µM).Cultures were supplemented with an antibiotic/ antifungal cocktail, which was reduced to 50% concentration on day 2 and discontinued on day 4 of culture.Culture medium was changed on a Monday-Wednesday-Friday schedule, and the plates were monitored for basal cell colony development.Cultures were discarded if colonies were not observed on culture day 14.Successful cultures (p0) were expanded at passage 1 (p1) and cryopreserved as frozen stocks at passage 2 (p2).

Clone-forming cell frequency assay
Clone-forming Cell Frequency (CFCF) was determined using a limiting dilution method that has been previously described [21] and modified for use with airway tissues [14].Passage one cells were distributed to wells of a 96-well plate using two-fold dilution from 100 cells/ well to 1 cell/well.Six replicates were seeded for each cell input.Following 9 days of culture, cells were fixed with NBF, and stained with Giemsa stain.Wells were scored as positive or negative for the development of colonies.Linear regression analysis was used to calculate the CFCF.

Population Doubling Level calculation
Population Doubling Level was calculated as: PDL = 3.32 (logXe -logXb), where Xe is the basal cell number obtained at passage and Xb is the basal cell number plated (https:// www.atcc.org/ resou rces/ cultu re-guides/ animal-cell-cultu re-guide).Cells were stained with trypan blue and counted using a hemacytometer [15].

Relative telomere length determination
Telomere length was determined based on previously published methods [22,23] with minor modifications [17].Specifically, cells were recoverd using the double trypsinization method.Genomic DNA (gDNA) was purified from 50,000 cells using the DNeasy Blood & Tissue Kit (Qiagen Cat# 69506) and quantified using a Nanodrop 2000 spectrophometer (ThermoFisher, Welmington, DE).Both the telomere and 36B4 single copy gene control assays used 40 ng gDNA/reaction.Targets were amplified using previously reported standards and primers [22].For the telomere assay, the magnesium concentration in the PowerSYBR Green Master Mix (Applied Biosystems Cat# 4367659) was lowered to a final concentration of 0.9625 mM through the addition of 0.5 M EDTA (0.0385 µL per 20 µL reaction).Previously published cycling parameters for telomeres [23] and 36B4 [22] were used.Jurkat (ATCC) and U1301 (human T-cell leukemia cell line (Sigma Cat# 01051619) were used as short and long telomere controls, respectively.A standard curve was generated for each experiment and was used to insure a linear relationship between DNA input and telomere length.Prior to analysis, the telomere standards were incubated at 95 °C with shaking for 5 min, and immediately placed on ice.Three technical replicates were generated for each standard curve point and biological sample.Relative telomere length was calculated according to the Telomere/Single copy gene control (T/S) ratio.

Correlation of clinical and culture phenotype
Univariate chi-squared and t-test analyses were used to compare patient demographic and clinical characteristics between those who did and did not have squamous basal cells.Clinical characteristics included prematurity; significant pulmonary, cardiac, gastrointestinal, and/or neurological history; known chromosomal abnormalities; congenital defects in the head and neck (including supraglottic, glottic, subglottic, tracheal stenosis); and history of airway surgery.Prematurity was coded as an ordinal variable per World Health Organization guidelines, with a moderate to late preterm group (gestational age between 32 and 37 weeks), a very preterm group (28 to 32 weeks), and an extremely preterm group (less than 28 weeks).

Statistical analysis
All statistical calculations were conducted with Stata 17.0 (College Station, TX).Summary data are presented as the mean and standard deviation.Symbols represent the values from individual donors.Sample size is noted in the figure legends.Normally-distributed data were assessed by student's t-test or 2-way ANOVA.Skewed data sets were assessed using the ROUT method and analyzed using GraphPad Prism, version 9.0.0 (GraphPad Software Inc., CA).A significance threshold of 0.05 was used.Specific analysis details for each data set and corresponding N for each experimental set are indicated in the figure legends.

Power analysis
Power analyses were completed with our pediatric population (total sample size of 48 patients), demonstrating that with our univariate chi-squared analyses with an error probability of 0.05 and power of 0.8, we would be able to detect a large effect size (r = 0.49).In our adult population (total sample size of 6 patients) with similar power analyses, we would be able to detect an effect size of 0.78.In order to detect a medium effect size (r = 0.3), we would need a total sample size of 143 patients in each of the pediatric and adult cohorts.

Adult donor characteristics
Six adult donors were recruited for this study.The cohort was 50% (N = 3) male, with mean age of 55 years (SE 13.2 years).50% (N = 3) of the patients underwent direct laryngoscopy (DL) for subglottic stenosis and 33% (N = 2) for tracheal stenosis.The final patient had a normal airway and was recruited as a healthy control.Additional patient characteristics of the adult population are detailed in Table 2.

Basal cell number and phenotype
Cytospin slides were generated from uncultured cells and stained for a pan-basal cell marker keratin 5 (K5) and a marker of basal cell activation keratin 14 (K14) (Fig. 1A).Regardless of age (pediatric vs. adult) or region sampled (normal vs. stenosis), basal cell number (K5+) and activation status (K5/K14+) were similar (Fig. 1B).These data indicate that the brush biopsy protocol recovered similar numbers of basal cells and that these metrics did not vary by age or region (Fig. 1C).
As the brush biopsy contains a heterogenous population of epithelial cells including differentiated mucousproducing and ciliated cell types, cell isolates were cultured using the mCRC method which selects for basal cells.To determine if age or the region of airway stenosis sampling influenced basal cell expansion in vitro, cytospins were generated from passage 2 cells and stained for K5 and K14 (Fig. 1A).The percentage of basal cells increased with passage.Basal cell number (K5 + only) and activation state (K5/K14+) did not vary as a function of age (Fig. 1B) or airway region sampled (Fig. 1C).

Regenerative capacity
To estimate the therapeutic potential of cells that are recovered by brush biopsy, we determined the clone forming cell frequency at passage 1.This study used the limiting dilution assay and results were reported as the Clone-forming Cell Frequency (CFCFx1000, Fig. 2A).The mean CFCF was 195 +/-210 (n = 33 donors).This  value was lower than the CFCFx1000 for basal cells that were isolated by digestion of explanted tracheal tissue [24].The pediatric and adult groups and the normal or stenotic sampling sites had similar numbers of clone forming cells.
To determine if these clone forming cells could be amplified, they were passaged and cell yield was determined (Fig. 2B).These results were benchmarked against a scaffold that would be used for the reconstruction of a pediatric airway [25].At passage 1, all samples contained sufficient cells to cellularize a 3.5 by 0.5 cm 2 graft scaffold at low cell density (~ 7000 cells/ cm^2) and 98% could cellularize a scaffold at high cell density (~ 100,000 cells/cm^2).At passage 2, all samples contained sufficient cells for both recellularization models.
To determine if these cells had reserve regenerative potential, we first calculated the number of populations doublings [15].At passage 1, most samples underwent 0-5 population doublings (Fig. 2C).However, cell number did not increase in 13 samples.At passage 2, all samples underwent 4 population doublings.The total number of population doublings was ~ 25% of basal cell life-span [20].
Finally, we used telomere length as a surrogate for biological age [26].DNA was recovered from cryopreserved p2 cells (48 pediatric and 5 adult subjects).Telomere length was analyzed using a modified PCR assay and is presented as the relative telomere length [17].Since telomere length varies as a function of chronological age and has been associated with lung diseases [27], we evaluated pediatric and adult samples separately.At passage 2, the relative telomere length for pediatric samples was 2.3+/-0.42,and 41% of samples were in the short telomere range (the samples with a value </=1, Fig. 2D).Relative telomere length in adult samples was 1.6+/-0.3and most samples were in the short telomere range (Fig. 2D).No significant differences were detected between normal or stenotic samples (Fig. 2E).Collectively, these data indicate that the airways of pediatric and adult 54 donors have sufficient regenerative capacity to recellularize an airway scaffold.

Alteration of basal cell identity: squamous basal cells
Culture morphology was qualitatively assessed at p1.Typical basal cells (normal phenotype) formed discrete colonies on the fibroblast feeder layer that were generally spherical or ovoid in shape and had distinct edges.Normal phenotype colonies required extensive trypsin exposure to remove them from the culture plate and to dissociate the cells into a single-cell suspension for passage and handling.Most pediatric donor collections yielded basal cell colonies of typical phenotype (59 of 69 samples, 85.5%, Fig. 3A).These cells expressed high levels of K5 and K14 (Fig. 3B, C).
In contrast, squamous basal cells (Fig. 3F) were observed as morphological variants at p1.These colonies exhibited a diffuse lattice-like structure with poorly defined boundaries and expressed lower levels of K5 than typical cells (Fig. 3H, I).Both the number of K5 positive cells and the intensity of the staining were decreased relative to cells recovered from typical morphology colonies.Abnormal cells grew more rapidly than typical basal cells and were easily dislodged with trypsin.Ten cultures originating from seven pediatric patients (14.5%) contained squamous basal cell colonies at passage 1.Samples containing these abnormal cells were excluded from the previous data sets (Figs. 1, 2 and 3).
Cell phenotype was further evaluated using flow cytometry.Cells with a typical phenotype demonstrated high expression of the markers CD151 and CD49f (Fig. 3D, E), which are commonly expressed on epithelial basal cells [13].In contrast, cells with the squamous basal cell phenotype had decreased or absent expression of these markers (Fig. 3K-J).Squamous basal cells from pediatric and adult donors exhibited similar morphology and marker expression profiles.

Table 3 Squamous basal cell phenotype associations with the pediatric population
Univariate analyses were conducted, and statistically significant predictors were calculated with p<0.05 and marked (**) a DLB: direct laryngoscopy and bronchoscopy Subglottic stenosis as assessed by direct laryngoscopy and bronchoscopy approached significance (50% vs. 17%, p = 0.06).
Squamous basal cells were observed in one adult patient with subglottic stenosis and two patients with tracheal stenosis.Adult samples containing these abnormal cells were excluded from the previous data sets (Figs. 1, 2 and 3).The presence of squamous basal cells did not correlate with direct laryngoscopy findings (subglottic stenosis 50% vs. tracheal stenosis 100% vs. control 0%; p = 0.2) including aspiration, obstructive sleep apnea (OSA), past airway surgery, or current/previous tracheostomy (all p > 0.05).

Discussion
Our study revealed that the airways of pediatric and adult donors have sufficient regenerative capacity to recellularize an airway scaffold.However, we identified squamous basal cells in a subset of donors and associated their presence with chronic lung disease and tracheostomy in the pediatric population.
The basal cell is the tissue-specific stem cell of the conducting airway epithelium [28].Common pediatric airway diseases including bronchopulmonary dysplasia [29] and laryngotracheal stenosis [30] can result in epithelial injury from mechanical ventilation, direct trauma, infection, and inflammation.Repeated injury causes selective activation of basal cells [17] and promotes proliferation.These events can lead to terminal differentiation and squamous metaplasia [31][32][33].
In conditions such as chronic obstructive pulmonary disease (COPD), pathologic remodeling of the airway epithelium includes squamous metaplasia.This condition is likely to be the result of environmental exposures and a chronic inflammatory response and is thought to contribute to poor ciliary function and mucus clearance [28].Sustained pathogeneic squamous phenotype has been associated with the pathophysiology of tracheobronchial disorders with unknwon etiology, accompanied by recurrent infections, nodules formation, and absence of normal ciliated respiratory epithelium [34].Squamous metaplasia has also been reported in the proximal airways following intubation and tracheostomy placement.Although the clinical implications of squamous metaplasia remain unclear [13], McKeon's group identified a squamous basal cell subtype that self-renewed and retained its phenotype across multiple passages and after transplantation into immunocompromised mice [31].We identified similar cells in our cohort and report that they are more prevalent in pediatric patients with bronchopulmonary dysplasia, subglottic stenosis, and those who underwent tracheostomy placement.
The association between bronchopulmonary dysplasia and squamous basal cells requires further study.The observed correlations may be a consequence of direct trauma to the airway during intubation and surgery in combination with free radical changes from prolonged intubation.It is also possible that the squamous metaplasia phenotype is congenitally predisposed in certain pediatric patients, although the relationship was not clearly demonstrated [13].If substantiated, this predisposition would explain why premature patients have significantly higher rates of squamous basal cells compared to their full-term peers [35].

Future directions
The rise of tissue engineering has provided promising techniques with the intention of creating viable tracheal replacement grafts [8,36,37].Successful grafting not only requires an optimal cell source and scaffold for integration but also the appropriate cellular microenvironment for complete epithelialization [8].Although our data indicate that brush biopsy cells have significant therapeutic potential, a subset of samples contained squamous basal cells that do not regenerate the epithelium.Further studies could determine if cell rejuvenation methods would improve the performance of samples with low-performance metrics.

Study limitations
We recognize that our understanding of the squamous basal cell phenotype is limited by sample size and that the study may be underpowered to detect smaller effect sizes, particularly those related to tracheal stenosis.We also acknowledge that our cross-sectional study design could be enhanced by a longitudinal study.Finally, we understand that our sampling method may have underestimated the frequency of squamous basal cells in pediatric and adult airways.

Fig. 2 Fig. 3
Fig. 2 Regenerative capacity of pediatric and adult basal cells.A Bar graph represents the Clone-forming Cell Frequency of pediatric and adult basal cells after 1 passage in culture (n = 30).B Prediction of the number of cells recovered after 1 or 2 passages in culture (n = 48) able to cellularize a scaffold at low density (LowD) or high density (HighD).C Bar graph of the in vitro population doubling after 1 or 2 passages in culture (p1 n = 35, p2 n = 35).D Relative telomere length of the basal cells of pediatric (n = 48), and adult (n = 6) donors divided by E Normal (n = 59) and stenotic (n = 18) region.Bars represent mean ± S.D. Statistical analysis was performed using paired, 2-tailed Student's t-test

Table 1
Clinical characteristics of pediatric patients

Table 2
Clinical characteristics of adult patients